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Abstract. The paper is concerned with the maximum likelihood estimator (MLE) 
of the unknown drift parameter 9 G R in the continuous-time regression model 

X t = 6t + B t +Bf, te[0,T] 

where B t is a Brownian motion and Bf is an independent fractional Brownian 
motion with Hurst parameter H € (|, 1). We derive the exact formula for the 
MLE in terms of the solution of an integral equation and find the asymptotic 
distribution of the estimation error. In particular, it turns out that the Brownian 
part does not contribute to the asymptotic variance of the MLE. 

Another contribution of this paper is a formula for the Radon-Nikodym de- 
rivative of the probability, induced by the mixed fractional Brownian motion 
£t = Bt + Bf , H > 3/4 with respect to the Wiener measure. 



1. Introduction and the main result 

Consider the continuous-time regression model 

X t = 6t + aB t + Bf, te[0,T], (1.1) 

where B t is a Brownian motion and Bf is an independent fractional Brownian 
motion (fBm) with Hurst parameter H E (|,l), i-e., zero mean Gaussian process 
with the covariance function 

EBfBf = U\ t \ 2H + \ s \ 2H -\t-s\ 2H y s,te[0,T}. 
As is well known, for He (|, 1) the process Bf exhibits the long-range dependence 



property 



oo 

£ 

3=1 



EBf(Bf +1 -B?)=oc, 



and hence £ t := aB t + Bf, called in [1] the mixed fractional Brownian motion (fBm), 
can be thought of as observation noise with both "white" and heavily correlated 
components. The mixed fBm has a number of peculiar probabilistic properties, 
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studied in e.g. [I], [2], [T7], which have some relevance to mathematical finance (see 
e.g. 0). 

The constant a > 0, controlling the intensity of the Brownian part, and the Hurst 
parameter H can be reconstructed precisely from the trajectory X T := {X t ,t G 
[0,T]} (see, e.g., [1]) and hence are assumed to be known. 
Given the sample path X T , it is required to estimate the unknown drift parameter 
£ K. The parameter estimation problems in models with mixed fBm have been 
considered in the recent monographs [5] and [T2], where the construction of the 
maximum likelihood estimator (MLE) of 9 appears as an open problem (see Remark 
(iii) page 181 in [TJ] and the discussion on page 354 in [H]). Our main result aims 
at filling this gap: 

Theorem 1.1. The MLE of 9 is given by 

£ 3 ( t ,T )d X 
/„«(«, T)dt 

wher^ the function g(t,T), t G [0,T] is the unique L 2 [0,T] solution of the integral 
equation 

r T 

a 2 g(t,T) + H(2H- 1) / g(s, T)\s - t\ 2H ~ 2 ds = a, for a.a. t G [0,T]. (1.3) 



This estimator is strongly consistent and the corresponding estimation error is nor- 
mal 

§t — 9 ^ N ( , — if, — — ), (1.4) 

y^g{t,T)dt)' 

with the asymptotic behavior of the variance 

9 9 ff r „ 2HT(H + l)r(3 - 2H) , 
hm T 2 ~ 2H E e (9 T - 9) 2 = X H := K t{ 3_ _h) ' ( L5 ) 

where T(x) is the standard Gamma function. 
Remark 1.2. 

(1) The asymptotic variance in ( II. 5p is independent of a and coincides with the 
asymptotic variance of the MLE in the problem with a = 0, i.e., estimating 
the drift of fBm without additional Brownian component (see Section 5.1 in 
[B]). This means that the Brownian part is asymptotically negligible. 

(2) Actually, it will appear (see Remark 12.31 below) that for some constant kh 
the following limits hold: 

l im = fcjf(T) t) k - h H\- h (T - t)^ H , (1.6) 



stochastic integral in the numerator is defined through the usual limit procedure, recalled 
in Subsection! 
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The MLE weights functions: mixed fBm and fBm 



3 
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Figure 1. The MLE weight function for mixed fBM versus fBm 
(a = 1, T = 1, H = 3/4) 



lim / g -^Hdt = [ k H (T,t)dt = \ H 1 T 2 ~ 2H . (1.7) 
Jo cr J 

Hence, the limiting form of ( 11.21) for o tending to is 

Or = j^h j\ H (T,t)dX t , (1.8) 

which is nothing else but the expression obtained in [6] for the MLE in the 
model with o = 0. It is easy to check that this estimator is applicable to 
the data X, generated by the model with any o > and its asymptotic 
variance coincides with (jl.5p . In other words, the estimator ( 1 1 . 8 [) has the 
same asymptotic accuracy as the genuine MLE. 



Remark 1.3. 

(1) The proof of Theorem 11.11 suggest s an approximation procedure for the func- 
tion g(t, T) (see ( 12. 5ft and ( 12. 6p ). Its typical form, depicted in Figure [TJ versus 
the weight function from the estimator ( II. 8p . indicates significant difference 
in the non-asymptotic regime. 

(2) The integral equation ( II. 3p is known as the second type Fredholm equation 
with weakly singular kernel (see |11] ) or as the Wiener- Hopf equation on 
the finite interval. Its solution can be reduced to a particular instance of 
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the Riemann boundary value problem, which unfortunately doesn't seem to 
be helpful in our case. It is well known, however, that (ll.3p has a unique 
continuous solution, which enjoys some regularity properties (see, e.g., [T6]). 



Another interesting outcome of our approach is a formula for the Radon-Nikodym 
derivative of the probability, induced by the mixed fractional Brownian motion with 
respect to the Wiener measure. It is shown in [4] (see also [2]) that for H > 3/4 
these measures are mutually absolutely continuous, however no expression for the 
corresponding derivative is given. The derivative is calculated in terms of reproduc- 
ing kernels in [17], but the author points out that it might be hard to obtain more 
explicit expression (see remark (iii) on page 63). The following proposition gives a 
representation formula in terms of the solution of (11.31) : 

Proposition 1.4. For H > 3/4, the probability ^fi, induced by £t := B t + Bf , 

t £ [0, T] is absolutely continuous with respect to the Wiener measure [i w and 



dfi 



with tp a (£) = f° fjj^r and g(r,t) = §- t g(r,t) where g = (g(s,t), < s < t < T) 
satisfies the equation 

g(s,t) + H(2H -I) [ g{r,t)\r - s\ 2H ~ 2 dr = 1, s £ [0, t\. (1.9) 
Jo 

2. Proof of Theorem 11.11 

2.1. The likelihood function and the MLE. Let B = (B t ) t£[0tT] and B = 
(i? t // ) te [o i T] be processes defined on a measurable space (fl, 5F) and Fg be a probabil- 
ity, under which B and B H are independent, B H is a fractional Brownian motion 
with Hurst parameter HE (|, 1) and B is a Brownian motion with drift -, i.e., 

aB t = 9t + aB t , te[Q,T\. 

Under Fg, the process X = oB + B H is the mixed fBm with drift 9 as defined in 
(11. ip . By Girsanov's theorem and independence of B and B H 

d¥ e (b~ ie 2 

dF =eXP [a Bt ~2^ T 

In order to eliminate some confusion, existing in the literature (see, e.g., [Hj), let us 
stress that this derivative is not the likelihood for the problem at hand, since ,e.g., 
it is not measurable with respect to the observed cr-algebra 3? = a{X t ,t £ [0, T]}. 
The proper likelihood function is obtained by conditioning on 9^ . 

More precisely, let \xq be the probability induced by X on the space of continuous 
functions with the usual supremum topology under probability ¥$. Then for a 
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measurable set A, 

dW / dW 3 \ 

fie(A) =F S (X eA)= E ^l {XeA} = E E (^|9f ) MxeA} 

$(x)/i (<ix), 



where is a measurable functional, such that 



$(X)=E (^|3f), P -a. S . 



The latter means that /!# <C fiQ for any S 6 1 and, since B = B under P , the 
corresponding likelihood function is given by 

W) =H,(^|3#) = E„ (exp ( 6 -B t -^t) = 

exp (V + I^(^-T)). 

The latter equality holds with M T := E (B T \^) and Vt = E (B T - M T ) 2 , since 
the process (B,X) is Gaussian and hence the conditional distribution of Bt given 
9^ is Gaussian as well. 

Let (Ht) and (S^) be the natural filtrations of (B, B H ) and X respectively and 
set 

M t = Eo(B t |3f), te[0,T]. (2.1) 

Since B is an (^FJ-martingale and 9y C 3^, the process M is an (3 r ^")-martingale 
with respect to P . Moreover, since V t = E (-B t 2 |3 r f ) - M t 2 and B} - t is an (3^)- 
martingale, for s < t, 

E (M 2 -{t- V t )\<r?) = E (E (B?\J?) - t\3*) = E (5 2 - f|3* ) = 
E (B 2 s \?*)-s = M*-(s-V s ), 

i.e., the quadratic variation process of the martingale M is (M) t — t — Vt, and the 
likelihood function reads 

L T (X; 0) = exp f-M T - ~(M) T ) . 
The MLE of being the maximizer of the above expression, is given by 
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This estimator is unbiased: 

E e a— f- = aE L T (X; n 



(M) T u 1K ' ' (M) 



T 



ff (M)7 eXP ("2^ (M>r )^ eXP (2^ <M>r|=( ' 
with the variance 

, 2 



E e (e T -ey = E e e 2 T -e 2 = <t 2 e*^i -e 2 = 

^e X p(4^(M) T ) Eo e X p(V)^-^ 
exp f-l^(M) T ) -^E exp ( -M T ) - 6 



* 4 f«W + <^_'- ff2 



(2.2) 



•"> 



To recap, the MLE error is a zero mean Gaussian random variable with variance 
a 2 / (M)t- Next we shall derive an explicit characterization of the martingale M 
in terms of the solution of the integral equation (jl.3p and will find the appropriate 
asymptotic as T — » oo. 

2.2. The martingale representation. Let us recall briefly some relevant prop- 
erties of the integrals with respect to fractional Brownian motion. Following the 
notations of [10], define the spaces 

L 2 [0,T] := {/ : [0,T] ^ R such that J f 2 {u)du < oo}, 

|A|J~s : = j/ ; [o,T] ^ R such that jf jf |/(u)||/(t;)||u - v\ 2H ~ H dudv < oo}, 

A^H := |/ : [0,T] R such that ^ [s^' H (l^u H ~^ f{u)){s)^ ds < oo}, 
i 

where 1 T _ 2 is the Riemann-Liouville fractional integral operator (see [IS])- For 

if G (|, 1) the inclusions L 2 [0, T] C |A|J^ C A^ hold (see Remark 4.2 in [ID]). 
For the simple function of the form, 

n 

f( u ) = 5^ /fc G M > = ni < n 2 < ... < n fc = T, 

fc=i 
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the stochastic integral with respect to B H is defined by 

n 



f 
Jo 



k=l 



H-- 

Since the simple functions are dense in A T 2 (see Theorem 4.1 in [10]) , the definition 
of j Q f(t)dBt is extended to / £ A T 2 through the limit 

f f{t)dB? :=lim f f n (t)dBf, 
Jo n Jo 

where /„ is any sequence of simple functions, such that lim n ||/ — / n || H _i = 0. 



It turns out however (see Section 5 of [H]]), that the image of A r 2 under the 
map / !->• f Q f{t)dB^ is a strict subset of sp^ T ^B H ), the closure in L 2 (r2,3 r , P ) 
of all possible linear combinations of the increments of B . In other words, some 
linear functionals of B H cannot be realized as stochastic integrals of the above type. 
Nevertheless we have the following: 

Lemma 2.1. Assume H £ (|, 1) and let r\ be a Gaussian random variable, such 
that (r],X t ), t £ [0,T] is a Gaussian random process. Then there exists a function 
g(-,T) £ L 2 [0,T], such that 

EoM3f) =E r/+ / g(s,T)dX s , ¥ -a.s. (2.3) 
Jo 

Proof. Following the arguments of the proof of Lemma 10.1 in [TJ, let (tj), i = 0, 2™ 
be the dyadic partition of [0, T], i.e., U = i1~ n , i = 0,...,2 n and 3 r ^ n = o~{X ti — 
X u ^,i = 1, 2 n }. Then 3^ n / j£ and by the martingale convergence 

UmEoM^J =Eo(77|ff£), P - a.s. (2.4) 

n 

as well as in L 2 (f2, 3 r , Pq), since Eo(?7|3 r ^ n ) are uniformly integrable. By the Normal 
Correlation theorem, 

2" 

Eofo|3*J =E 77 + f>?_ 1 (X, i -X ti _ 1 ), (2-5) 



i=l 



with constants i = 1, 2 n . Define 



2" 

9n.it, T) := ^2 9i-il{[U-i,ti)}(t), (2.6) 

i=l 

then 

Eo(f7|2£ n )=Eo77 + <7 / g n (t,T)dB t + [ g n {t,T)dB? , 

Jo Jo 
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and 

E 



2 P 

o (e (77 1 5F^ n ) - Eo (77 1 9^ m ) ) = a 2 ^ (g n (t,T)-g m (t,T)) 2 dt+ 

°h / (g n {t,T) - g m (t,T))(g n (s,T) - g m (s,T))\s -t\ 2H ~ 2 dsdt, 
Jo Jo 

where ch '■= H(2H — 1). Since the kernel in the last integral is positive definite 
limsupa 2 / (( 7n (t,T)-( ?fn (t,r)) 2 ^<limsupEo(E (r ? |^ n )-E (r ? |^ m )) = 0, 

n m>n Jq n m>n ^ ' 

where the latter equality holds by (\2A\i , since L 2 (Q, 5F, P ) is complete. Since L 2 [0, T] 
is a complete space, there exists a function g(t,T) G L 2 [0,T], such that lim n \\g — 

g n h = 0- Tnen 

E (E (r)\^)-E Q r)-a g{t,T)dB t - J g{t,T)dB^ < 

2 P*~^ 

3E (Eo(r / |^)-E (r / |^j) + 3a 2 J (g n (t,T) - g(t,T)) 2 dt+ 

3c H / / (<?„(*, T) - g(t,T))(g n (s,T) - g(s,T))\s -t\ 2H ~ 2 dsdt 0, 
Jo Jo 

H-- 

where the latter convergence holds, since L 2 [0,T] C |A| T 2 . □ 

Applying the above lemma, we obtain the claimed formulas ( jl.2[) and (II ,4ft : 

Lemma 2.2. Zet (M t ) be the (Uf) — martingale defined by (12. ip . T7ie following 
representations hold: 

M T = [ T g(t,T)dX t , (2.7) 
Jo 

and 

(M)t = o [ g(t,T)dt, (2.8) 
Jo 

where g(t,T) is the unique solution of (11.31) . 

Proof. By Lemma [2. 1| there exists g(-,T) G L 2 [0,T], such that 

M T = Eo(Br|2r) = / g(t,T)dX u F -a.s. 
Jo 
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holds. For an arbitrary function h G L 2 [0,T], 
E (b t -J^ g{r,T)dX^j jf h(s)dX s = 

E Qf dBt-a g(t, T)dB t - g(t, T)dB* \ [a h(t)dB t + jf h{t)dBf\ 
£ h{s) (a - a 2 g(s, T) - c H £ g{r, T)\s - r\ 2H ~ 2 dr^ ds. 

By the orthogonality property of the conditional expectation and by arbitrariness 
of h, it follows that g(t,T) satisfies (11.31) for almost all t G [0,T]. This solution is 
unique (see, e.g., [IS])- Further, since M is a Gaussian martingale, 

r-T \ 2 



(M) T = E M| =®o(J 9(8, T)dX s 
^ g(t,T)(a 2 g(t,T) + c H g(s,T)\s - t\ 2H ~ 2 ds^dt = a jf ^(t,T)dt. 



□ 



2.3. The large sample asymptotic. Finally we shall derive the asymptotic an- 
nounced in fO]) . Let /i := T 2H ~ X and define ^(u) := T 2H ~ x g{uT,T), u G [0,1]. 
Then (II. 3p reads 



/i 

and, moreover 



1 r 1 

AM + c h / 0»|u-Tf H " 2 dv = <r, we [0,1], (2.9) 
M Jo 



(M) T = a [ g(s,T)ds = aT 2 ~ 2H [ g^{u)du. (2.10) 
Jo Jo 

Define the operator K 

Kf{u) = c H f f(v)\u - v\ 2H ~ 2 dv, f G |A|£~* (2.11) 
Jo 



and the scalar products 

(/, h) :-- 

and 



f f(s)h(s)ds, f,heL 2 [0,l} 
Jo 



(f,h) K :=c H f [ f(v)h(u)\u-v\ 2H - 2 dvdu, h,fe\A\" 
Jo Jo 

In terms of these notations, the equation (12. 9p becomes 

9n + Kg^ = a. 
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We shall also consider the first type auxiliary integral equation 

Kg = <x, (2.12) 
which admits the unique solution (see Lemma 3 in [6] and the references therein): 

where P(a,i) = r r ( ff ■ 

The function 5^ := — g satisfies 

—b^ + = g. 

H-- 

Since g G £ 2 [0, 1] C |A| : 2 , multiplying by 5^ and integrating we obtain 

-\\S,\\l + \\SJ 2 K = -\(g,5,}\, 



and, in particular, U^IH < \{g, 5^) . On the other hand, by the Cauchy-Schwarz 
inequality |(<7, (S^)! < H^l^ll^lh and hence H^Ha < IMh- Note that <5 M also satisfies 



a 2 



g„ + KSp = 0. 

Multiplying both sides of this equation by g and integrating, we get 

2 

— (g ll ,g) + (K6 li ,g) = 0. 

But 

|2 ^ II r II II „|| , II „||2 ^ oil 1 1 2 



< \{5^g)\ + \\g\\i < \\5J 2 \\9h + Ml < 2\\g\\i < °o 



and hence 



= \(5„Kg)\ = \{KS„g)\ = -\(g„g)\ < j2\\g\\l ^ 0, 



and 



lim / gJu)du = / g{u)du. 
Jo Jo 

Finally, by the formulas (E2} and f l2~T0l) 

T 2 " 2H E e (^ - fl) 2 = T 2 " 2H = ^ ^> — ° 



(M) 
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The asymptotic ( 11. 5ft now follows, since 



9{r)dr= C 2 (H)T(2-2H) l -0*-** 
1 0(2-2H,H-l) 1 o ,3 TT 3 



-B{--H,--H) 



H(2H-1) T 2 (H-\) Y{2-2H)' y 2 '2 

i r(2 - 2#)r(ff - 1) i r 2 (|-g) 

H{2H - 1) T(f - H)Y 2 (H - I) T(2 - 2#) T(3 - 2#) 



2 ' V- 1 - 1 2' 

r(§-#) r(|-£0 



# (2iJ - i)r(# - |)r(3 - 2if ) 2#r(# + |)r(3 - 2H) ' 

where we used the property T(x + 1) = xr(x), x > 0. 

Of course, the martingale property of the process '■= g(t,T)d^ t , T > 
with £t '■= crBt + Bf and the representation of the error: 

M? 



6t — 



T 



(M H ) T > 

where {M h )t —> oo when T — ?• oo implies the strong consistency of 9t due to the 
strong law of large numbers for martingales. 

Remark 2.3. By means of a study similar to the previous, one can check that the 
limits dLSD and flTT]) hold. 

Remark 2.4. It is worth mentioning that, following the terminology of [9], [6] and [5], 
the martingale M H merits to be called the fundamental martingale associated to the 
mixed fractional Brownian motion £. It will play also a key role in the statistical 
analysis of models more general than (II. ip such as the mixed fractional Ornstein- 
Uhlenbeck process. The progress in this direction will be reported elsewhere. 

3. Proof of Proposition 11.41 

In the following lemmas, we shall first prove a number of useful properties of the 
solution of (11.91) . assuming hereafter that H > 3/4. 

Lemma 3.1. The solution g{s,t) of (11.91) is continuously differentiable at t > for 
any s E (0,t). The derivative g(s,t) := -j^g(s,t) satisfies the equation 

g(s,t) + c H [ g{r,t)\r-s\ 2H - 2 dr = -c H g(t,t)\t-s\ 2H - 2 , se(0,t), t > 0, (3.1) 
Jo 

and g 2 (s,t)ds < oo. 

Proof. The function gt{u) := g(ut,t), u G [0, 1], t > satisfies the integral equation 

g t {u) + cut 211 - 1 [ g t (v)\v - u\ 2H ~ 2 dv = 1, u G [0, 1]. 
Jo 
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This equation has a unique continuous solution for any t > (see [IE]) and, in 
terminology of [T3], any point A := t 2H ~ l is regular. Since for H > 3/4 the kernel 
(12.111) belongs to L 2 [0,1], it follows from ,e.g., Theorem on page 154 in [13], that 
the solution g t {u) is analytic at t > 0. By [16] the solution g t {u) is continuously 
differentiable at u G (0,1) and hence the function g(s,t) = gt(s/t) is continuously 
differentiable at t > for any s G (0,t). The equation (I3.ip . obtained by taking the 
derivative of both sides of (II. 9p . has a unique solution in L2[0,t] for H > 3/4. □ 

Lemma 3.2. The solution g{s,t) of (II. 9p is suc/i that g(t,t) ^ /or all t > 0. 

Proof. Letting g'(s,t) := -j^g(s,t) and taking the derivative of (II. 9p . we obtain 

g\s,t) = -c H ^- g{r,t)\r - s\ 2H - 2 alr = -c H ^- (J g{u + s,t)\u\ 2H - 2 du^ = 

/t—s 
g'{u + s, t) |n| 2f/ - 2 du + c H g(t, t) \t - s\ 2H ~ 2 - c H g(0, t)\s\ 2H - 2 = 

- c H £ g'(r, t) \r - s\ 2H ^ 2 dr + c H g(t, t)(\t- s\ 2H ~ 2 - \s\ 2H ~ 2 ) , 

where we used the obvious symmetry g(t — s,t) = g(s,t) and g(t,t) = g(0,t) in 
particular. Now suppose g{t,t) = for some t > 0. Then 

g'{s,t) + c H [ g'{r,t)\r-s\ 2H - 2 dr = 0, s G [0, t\. 
Jo 

This equation has the unique solution g'(s, t) = 0, i.e., g(s, t) is a constant function. 
But since g{t,t) = 0, it follows that g(s,t) = for all s G [0,t], which contradicts 

<ra. □ 

Lemma 3.3. The solution g{s,t) of (II. 9p zs suc/i that for all t > 0, 

/ g(s,t)ds= [ g 2 (s,s)ds. (3.2) 

Proof. We shall use Krein's method for solving integral equations with difference 
kernels on a finite interval (see §13.13-1 in [H]). Let y(s,t) satisfy the equation 

y(s, t)+ f y{r, t)k{r - s)dr = f{s), s G [0, t], 



where k{u) = ch\u\ 2H 2 and / is a continuous function and consider the auxiliary 
equation with parameter r G [0,t\: 



9(s,t)+ g(r,r)k(r - s)dr = 1, s G [0, r]. 
Jo 

Then ^ 

y(s, t) = t) - f g{s, r)F'(r)dr, (3.3) 
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where 

F (T) = ^T^T f 9(s,r)f(3)da. (3.4) 

In particular, by uniqueness of the solutions, for / = 1 we get y(s,t) = g(s,t), 
s G [0,t] and setting s := t in Q, we see that F{t) = 1. Now by (JOK 

# 2 (M) = -j t J g(s,t)ds, 

and (13. 2p follows by integration. □ 
Now we are ready to prove Proposition 11.41 By Lemma 12.21 and Lemma | 



/ g(s,t)ds = f g 2 {s,s)ds, t e [0, T\. 
Jo Jo 



(M)t 

Hence by the Levy theorem and Corollary [3] 

f l 1 

W t = / zdM s , t G 0, T], 

Jo 9{s,s) 

is a Brownian motion. On the other hand, 

M t = J g{s,t)d£ 8 = J g{s,s)d£ s + J (g(r,t) - 0(r,r))d£ r = 
/ g(s,s)d£ 8 + / / g(r,s)dsd£ r = / g(s, s)d£ s + / / g(r,s)d£ r ds, 

Jo Jo Jr Jo Jo Jo 

and hence 

W t = [ —t~ — rdM s = £ t + / / ip^-d&ds =: £ t + / <p,(Z)ds. 
Jo g{s,s) J J g{s,s) Jo 

The desired claim follows from Theorem 7.7 in [7], once we check 

/ E$(B)dt < oo and f E^ 2 (f)dt < oo. (3.5) 
Jo Jo 



Since <£>t(") is additive and £ t = B t + -B^, where i? and B H are independent, it is 

g(t,t) 



enough to check only the latter condition. By (13. ip the function R(s, t) • 



satisfies 



R{s,t) + c H [ R(r,t)\r - s\ 2H ~ 2 dr = -c H \t - s\ 2H ~ 2 , se(0,t),t>0. (3.6) 
Jo 
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Hence for H > 3/4, using the latter equation we get 

mp*(£) = E( f R(r,t)d£ r } = [ R 2 (s,t)ds + c H [ [ R(s,t)R(r,t)\r - s\ 2H ~ 2 drds 
\Jo J Jo Jo Jo 

jf R(s, t) (r(s, t) + c H 1^ R(r, t) \r - s\ 2H - 2 dr^j ds = -c H R(s, t) \t - s\ 2H ~ 2 ds < 

c H (/' RHs, t)d s) " 2 (/' | t - s^ ds ) = -0= (/' RHs, ,)*.) " 2 
Since the kernel is positive definite, multiplying ( 13. 6\\ by R(s, t) and integrating gives 

j\ 2 (s,t)ds<-c H j\(s,t)\t-s\ 2H - 2 ds< V 4 ^_ 3 ^R 2 (s,t)ds^j ' t 2H ~ 3 / 2 , 
and consequently 



t \ V2 

' 2 f e fU. <T ^ *2fT-3/2 



2 

Plugging this bound back gives E</?f(£) < 4 ^_ 3 t 4H ~ 3 and in turn 

r r 2 /-T 2 



E^(fWt < — ^— / t w - 6 dt = t ^ -T 

-4:H-3J (AH - 3) (AH - 2) 

which verifies (I3.5P and completes the proof. □ 
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